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1 Introduction 



In these notes we discuss the topological nature of some problems in condensed matter 
physics. This topic has been widely studied in various contexts. In statistical mechanics, 
the possible stable defects in an ordered system have been classified according to the nature 
of the order parameter (e.g. scalar, vector, matrix) and the space dimensionality of the 
system using homotopy groups 0. Then, the discovery of the quantum Hall effects and 
the role played by stable integers or rational numbers for systems with few or no conserved 
quantum symmetries have motivated several topological models of quantum condensed 
matter systems [0, |3], ffl . A combination of these two ideas of defects classification and 
microscopic quantum models has been used in the description of superfluid 3 He ||. 

Here, instead of trying an exhaustive review of problems where topological ideas may 
play a role, we present the basic constituents needed in a geometrical description and 
calculate the related topological numbers. The ideas and methods developed in mathe- 
matics and mathematical physics to solve problems in geometry and topology are pretty 
sophisticated and sometimes expressed in a way unfamiliar to the physicist. We adopt 
the language of differential geometry to present this subject, since it is adapted to de- 
velop some intuition towards more elaborate concepts like fiber bundles, connexions and 
topological invariants. 

In the last two sections, we shall discuss the problem of superconducting billiards 
within the Ginzburg-Landau approximation. This problem is interesting for several rea- 
sons. First, it is a non trivial example on which topological methods naturally apply to 
give an elegant solution to the calculation of the ground state energy. It is also a situ- 
ation for which one can address the question of transition between different values of a 
topological number controlled by the boundary of the system. This question is similar to 
the transition between quantum Hall plateaus. 

Such problems are not only of academic interest. Our motivation has been triggered by 
a set of new experimental results obtained on small size aluminium disks U [7|] in a regime 
where their radius R is comparable with both the coherence length £ and the London 
penetration depth A. The magnetization, as a function of the applied magnetic field, 
presents a series of jumps with an overall shape reminiscent of type-II superconductors, 
although a macroscopic sample of aluminium is a genuine type-I superconductor. These 
notes end with a theoretical analysis of these experimental results. 

2 Different iable Manifolds 
2.1 Manifolds 

A differentiable manifold M of dimension n is a space which looks locally like an open set 
of R n . On the vicinity of each point of M one can define a local coordinate system 
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where each point is represented by a set of n real numbers (x 1 , x n ). There exists geomet- 
rical spaces which are not differentiable manifolds like arithmetical ensembles (Q, Z/5Z), 
fractals, objects with branching points (Feynman diagrams). On the other hand vector 
spaces, spheres, projective spaces, matrix groups GL(n, R), SO(2), SO(3) (direct isome- 
tries), SU(2), SU(n) (complex isometries) are differentiable manifolds. In order to do 
differential calculus on a manifold, different coordinate systems are required to be com- 
patible with each other: the local transformations between one system of coordinates to 
another have to be smooth and invertible. A quantity attached to a manifold is said to 
be geometrical (or intrinsic) if it is independent of the choice of a coordinate system. 
Example: the sphere S 2 in Euclidian space R 3 can be endowed with local cartesian coor- 
dinates, or with spherical coordinates (latitude and longitude). None of these two systems 
is global, for instance the spherical coordinates are singular at the poles, although from 
a geometrical point of view all points on the sphere are equivalent. Therefore the spe- 
cial role played by the poles are an artifact of the spherical coordinate system. On the 
other hand the tangent plane passing through a point of the sphere is an intrinsic object, 
although its equation looks very different in cartesian and spherical coordinates. 

2.2 Differential forms and their integration 
2.2.1 Tangent space 

Consider a manifold M and p one of its points. Let 7 be a curve in M passing through p. 
In a local coordinate system, the curve 7 is given by 



For t — 0, 7 passes through the point p represented by the coordinates (x 1 (0), ...,x n (0)). 
A tangent vector to the manifold M at the point p is by definition a tangent vector 
to a curve passing through p. For instance, the curve 7 defines a tangent vector at p 



7: [-1,1] 



M 



t » (x\t),...,x n (t)) 



(1) 




t=0 



with components: 




(2) 



If one uses a different set of local coordinates (X 1 , ...,X n ) in the vicinity of p, one finds 
that the components of v are given by 




(3) 



The transformation law from one set of components to the other is 




(4) 
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We use Einstein's convention of summation upon repeated indexes. The transformation 
law of the components of v is inverse of that of the partial derivatives ^j. Hence in 'classi- 
cal' (XIX century) mathematical literature, a tangent vector was said to be contravariant 
and the transformation formula was used as a definition: a tangent vector is an object 
whose components transform according to (f|) under a change of the local coordinate 
system. 

The set of all the tangent vectors at p to all curves 7 included in M and passing 
through p is a vector space, called the the tangent space of M at p and denoted by 
T p M. It has the same dimension as the manifold itself. With the help of (U), one can 
define a tangent vector, as an invariant object that does not depend on the coordinate 
system: 

= (5) 



The vectors 




l,...,n (6) 



are a basis of the tangent space T P M. This basis is local and depends on the point 
p. Usually one does not write explicitly the dependence on p. However it is extremely 
important to keep in mind that two tangent spaces T P M and T p iM at two different points p 
and p' are different vector spaces. It is not possible a priori to compare (i.e add, subtract) 
two tangent vectors v and v' at p and p'. A simple way to understand that is to realise 
that the transformation law ([|) is point dependent. Therefore v and v' could by chance 
have the same components in a special coordinate system but different components in 
another system. So it would be a mistake to say that they are 'equal'. 

The total tangent space TM of M is the collection of all tangent vectors at any 
point of M. Hence TM = \J p T P M. There is a natural function it from TM to M called 
projection: to each v G TM is associated p = tt(v) which is the point of M at which v 
is tangent to M. It can be shown that TM is a manifold of dimension 2n: if v belongs 
to TM one needs n components to specify at which point p the vector v is tangent and n 
more components to specify v in T P M. 

By definition a vector-field is a smooth field of tangent vectors to the manifold M, 
i.e. to each point p is associated a vector F(p) tangent to M at p. In a local coordinate 
system a vector-field is expressed as 

F(x\...,x n ) = F%x\...,x n )^- (7) 

2.2.2 Forms 

A O-form is a smooth scalar valued function on the manifold M . A 1-form is a linear 
function on vectors. Consider first the case of R 3 . A vector v = (x, y, z) in R 3 is written 
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in a basis v = xi + yj + zk. Let wbea 1-form on R. By linearity, u(v) is given by 

ui(v) = xuj{i) + yuj(j) + zuj(k). (8) 

Thus the question of calculating oj(v) boils down to the question of evaluating it on the 
vectors of the basis k. In R 3 , we define the 1-form dx by dx(i) = l,dx(j) = and 
dx{k) = and equivalently for the 1-forms dy and dz. We have shown that 

uj = u(i)dx + uj(j)dy + u(k)dz . 

The triplet (dx, dy, dz) is a basis for 1- forms called the dual-basis of k). 

More generally, one can consider a field of 1-forms on R 3 i.e. a 1-form uo with coefficients 

that vary from point to point. It can be written as 

uj = A(x, y, z)dx + B(x, y, z)dy + C (x, y, z)dz . 

The action of to on a vector field F = f(x, y, z)i + g(x, y, z)j + h(x, y, z)k is 

u(F)(x,y,z)=Af + Bg + Ch. (9) 

A k-form is a smooth multilinear and antisymmetric function on k-tuples of tangent vec- 
tors to M, all of them tangent at the same point. This can be implemented by considering 
two vectors in the plane, v\ = x\i + yij and i>2 = X21 + • Let be a 2-form, then, 

<f>(vi,v 2 ) = xmHhi) + x iV2(f>(iJ) +x 2 yi<p(],i) + x 2 y 2 (f)(jj) (10) 
= (xiy 2 - x 2 yi)<j)(i,j) 

because the antisymmetry condition implies <p(i,i) = 4>{j,j) = and <p(i,j) = —<f)(j,i). 
The requirement that differential forms be antisymmetric comes from the fact that we 
need to keep track of orientation. 

Examples: 

1. In R 3 the 2-form dx A dy is defined by dx A dy{i,j) = l,dx A dy(j,k) = and 
dx A dy{k, i) = and equivalently for the 2-forms dy A dz and dz A dx. 

2. More generally, in R n , the dual basis of the canonical basis (ei, ...,e n ) is denoted by 
(dx 1 , dx n ). It satisfies dx % (ej) = Sj, where 5* is the Kronecker delta. One defines the 
fc-form dx 11 A dx 12 A ... A dx %k , such that if (v\,v 2 , ...,Vk) is a /c-tuple of vectors of R™, the 
quantity dx 11 A dx 12 A ... A dx lk (vi, v 2 , ...,Vk) is the k x k determinant of the components of 
the vectors Vi along the directions defined by (e^, e i2 ...e ik ). The geometric interpretation 
of this number is known: it is the volume of the projection of the parallelepiped generated 
by (v\,v 2 , ...,Vk) on the linear space spanned by (e^e^, ...,ei k ). 

From example 2, one can prove that the set of /c-forms (dx 11 A dx 12 A ... A dx lk ) with 
H < H < ■■■ < ik is a basis of the vector space of the /c-forms in R™. Hence the space 
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of fc-forms has dimension k ,^[ k y - In particular all n-forms on R n are proportional to 
(dx 1 A dx 2 A ... A dx n ), which is nothing but the determinant. 

On an general manifold M, one constructs fc-forms locally for each tangent space T P M. 
For instance, a 1-form uj on M is a smooth linear function on T p M for every point p 
on M. Similarly, a fc-form w on a manifold M is a smooth collection of fc-forms on each 
tangent space T P M. For each p one defines a local basis ((dx l ) p , (dx n ) p ) dual to the 
basis (ei = fgfr) , e n = f gf^) ) of T P M. Then any fc-form w on M can be written 
in a local system of coordinates as follows: 

w = A h ,...,i k dx h A dx 12 A ... A dx ik (11) 

where the are rea l functions. Once again, the p dependence is not explicitly 

stated. 



2.2.3 Wedge-Product 

A 2-form ip can be constructed by forming the wedge-product if) = UJ\ A U2 of two smooth 
1-forms uj\ and U2 via: 

ip(u, v) = ujx(u)uj2{v) - ujx(v)uj 2 {u) (12) 

for any vectors u and v. More generally, if is a fc-form and if) a Z-form one can construct 
a (& + Z)-form ^A^ which acts on (fc + /)-tuples of tangent vectors at a point, by antisym- 
metrizing correctly the product (f>Aif). We shall not write the explicit formula here. It the 
same as the one used to construct a fermionic (antisymmetric) wave-function of (k + /) 
variables starting with two fermionic wave functions of respectively k and / variables. One 
has in particular 

Aif) = (-l) kl if) A0 . (13) 



2.2.4 The exterior derivative 

The exterior derivative d generalizes the usual operations of vector calculus. The exterior 
derivative d transforms fc-forms into (k + l)-forms. Let for instance A be a 0-form (a 
function), its exterior derivative dA is the 1-form 

r) A 8 A 

dA = ^—dx l + ... + ^—dx 11 . (14) 
dx 1 dx n v ; 

If v is a tangent vector then 

f)A 

is the rate of variation of the function A in the direction of v. This quantity is usually 
denoted by (\7A).v or V$A. We have seen that a fc-form w can be written in a local basis 



w 



J2 A h ,...,i k dx h A dx i2 A ... A dx ik 
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where the Aj^...^ are real functions. The exterior derivative operates on w by acting on 
each of the coefficients Aj^...^ via For example, for a 1-form = Acfe + Bdy + Ccfc, 
one obtains 



d(f) = dA A dx + dB A dy + dC A dz 

= (C y - B z )dy Adz + (A z - C x )dz A dx + (B x - A y )dx A dy (15) 

An important property of the exterior derivative is that it gives when applied twice, 

d 2 <p = (16) 

for any /c-form 0. This follows from the Schwartz identity for partial derivatives: 

d 2 A d 2 A 
dxidxj dxjdxi 

There is a Leibniz rule for the exterior derivative of a fc-form <p and a Z-form ip 

d(<j)Aip) = {d(j)) A-0 + (-l) fc <M# (17) 

Examples: In the Euclidean space R 3 , a vector field F is defined at each point (x,y,z) 
by F = Ai + Bj + Ck, where A,B,C are smooth functions of the coordinates. And 
let / be a function. Its gradient is the vector field V/ = f x i + f y j + f z k. Similarly, 
the application of the exterior derivative d on the 0-form (function) / gives the 1-form 
df = f x dx + fydy + f z dz . Hence, to the 1-form df is associated the vector field V/. The 
action of the rotational on the vector field F gives another vector field 

V x F = (C y - B z )i + (A z - C x )j + (B x - A y )k , 

and the exterior derivative d operates on a 1-form <p = Adx + Bdy + Cdz to give the 
corresponding 2-form 

d(f) = (C y - B z )dy Adz + (A z - C x )dz Adx + (B x - A y )dx A dy . 

Finally, one can check that the divergence of a vector field corresponds to d acting on a 
2-form to generate a 3- form. In summary, gradient, curl and divergence result from the 
application of d to 0-forms, 1-forms and 2-forms respectively. The relations V x VA = 
and V.V x F = are simply a consequence of d 2 = 0. 

2.2.5 Closed and exact forms. 

If cf) is a differential form defined on a manifold M with the property d<fi = 0, then <p is 
said to be closed. If it has the property that = dip for some form ip on each point 
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in M, then is exact. It follows from these definitions and from d 2 = 0, that every 
exact form is closed. But the reciprocal is not true with the important exception that 
on a simply connected domain M, i.e. a domain in which every closed curve can be 
continuously deformed to a point through deformations that remain in M, every closed 
1-form is exact. We shall discuss later in more detail the 1-form 

u = V dx + X dy. (18) 
x z + y z x z + y 

It is closed (duj = 0) but not exact on M = R 2 \(0, 0). 



2.2.6 Integration of forms. 



A fc-dimensional submanifold (or fc-chain) of a manifold M of dimension n > k is a subset 
which can be parametrized with only k coordinates. In this terminology, a 0-chain is a 
point, a 1-chain is a curve and a 2-chain a surface, fc-forms are the right objects to be 
integrated over /c-chains. Indeed fc-forms were invented by Elie Cartan for this purpose! 
(Their similarity with determinants is not by chance: when changing variables in an 
oriented integral one must multiply the integrand by a determinant and fc-forms were 
built to produce such a determinant by a change of the local coordinate system.) We 
shall simply consider the integral of 1-forms over 1-chains in the Euclidean space R 3 . Let 
C be an oriented smooth curve (i.e. a 1-chain) parameterized by f(t) = (x(t),y(t), z(t)) 
for t £ I = [—1, 1] and let uj = Adx + Bdy + Cdz be a smooth 1-form on C. Then, 

to = [ (A(x, y, z)x\t) + B(x, y, z)y\t) + C{x, y, z)z\t)) dt (19) 
c J-i 

Let F = Ai + Bj + Ck be the vector field which corresponds to the 1-form to. We obtain 
the more familiar expression 

I u= I F.dr (20) 
Jc Jc 

which represents the work of F along the curve C. 



Exercise : The winding number 

1. Show that on the unit circle 5* 1 parameterized by r(8) = (cos 8, sin 6) in the plane, the 
1-form u given by fll8|) is such that f S i u — 2tt. 

2. Show that if 7 is a path connecting two points P and Q of the plane R 2 \(0,0) the 
integral / u mesures the difference of the polar angles of Q and P, the center of the polar 
coordinates being (0,0). 

3. Deduce from 2. that if 7 is a closed path that encircles n times the point (0,0), 
W(7, (0, 0)) = f uj = n. In particular if 7 is a closed path that does not encircle (0,0), 
this integral is equal to zero. 
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4. The mapping W(7, (0, 0)) defined from the space of closed curves in R 2 \(0, 0) to the 
set of rational integers Z is called the winding number. It allows to classify different 
type of curves in R 2 \(0, 0). It is a simple example of a topological invariant (see section 
2.3). Two curves with the same winding number are said to be homologous. 

2.2.7 Theorem of Stokes: 

Let M be a compact oriented smooth manifold of dimension n with boundary DM (pos- 
sibly empty) and let dM be given the induced orientation M, [|. For a (n — l)-form 0, we 



The integral J c uj of a k-form is said to be path-independent if the value of this integral 
depends only on the boundary dC of the oriented fc-chain C. This implies that J c uj = 
for every closed /c-chain. This property can be used to state the important result: a form 
uj defined on a manifold M is exact iff J c uj is path-independent on M. For instance, the 
1-form (|18|) is not exact. 

All these properties are generalizations of well-known results in vector calculus. For 
instance, a vector field F is conservative if F — Vf for some function / (the potential). 
Let uj be the 1-form associated to F. A conservative F corresponds to uj = df so that uj 
is exact. This implies that uj is closed and since duj corresponds to V x F, this implies 
V x F = as well. Electrostatics results from the fact that in R 3 a closed 1-form is exact, 
hence the electric field is conservative. Another consequence is the widely used result (e.g. 
in thermodynamics) that a 1-form = pdx + qdy cannot be exact unless || = ||. 

2.3 Topological invariants of a manifold: 
2.3.1 Motivations 

Two manifolds are said to be homeomorphic if there is a continuous mapping from each 
other, with continuous inverse mapping. A topological invariant is a intrinsic charac- 
teristic of a manifold conserved by homeomorphism. These invariants reveal important 
features and help to classify different types of manifolds. Topological invariants can be 
numbers, scalars, polynomials, differential forms or more general algebraic sets such as 
groups, or algebras. Their importance for condensed matter Physics was recognized in a 
seminal paper of Toulouse and Kleman [1] for the study of defects (vortices, nodal lines, 
textures anomalies) and their stability as a function of external parameters. The analysis 
of [1] depends on general characteristics of the system under study (e.g. dimensionality of 
space, nature and symmetries of the order parameter) and not on the precise form of the 
equations governing the system. Loosely speaking, different types of defects correspond 
to different (non homeomorphic) geometrical structures. Therefore topological invariants 
can help to distinguish between them. Of course this general scheme does not tell how to 



have 




(21) 
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compute the relevant invariant in a given problem. A nice example is the Aharonov-Bohm 
effect in an infinite plane where the relevant invariant is the winding number (defined in 
the exercise of section 2.2.6). An thorough study was done in [[H]]. In the last chapter 
of these notes we discuss two dimensional superconductors. There the winding number 
measures the circulation of the phase of the order parameter around the vortices. In the 
following paragraphs we merely give a taste of the vast subject of invariants in a manifold. 



2.3.2 The Euler-Poincare Characteristic 

The Euler-Poincare characteristic x(M) of a manifold M is the oldest and the most 
celebrated topological invariant. We shall explain how to calculate it for a surface. Any 
given surface S can be tiled by triangles (of different sizes and shape). Such a partition 
is called a triangulation. For any triangulation of S, denote by V, E and F the number 
of vertices, edges and faces respectively. Then x{&) * s defined by 

X (S) =V-E + F . (22) 

This number does not depend on the chosen triangulation and is a characteristic of the 



surface S 11 



Exercise: Show for the sphere S 2 that x{S 2 ) = 2. Show for the torus T 2 that x{T 2 ) = 0. 
For a three dimensional manifold one needs to use tetrahedra to perform a 'triangulation' 
and the formula for x becomes x = V — E + F — T, where T is the number of tetrahedra. 
It is possible to generalize this notion to higher dimensional manifolds [H. 



2.3.3 De Rham's Cohomology 

We have seen (section 2.2.5) that any exact form is closed but that the reverse is not true. 
Forms which are closed but not exact reveal important topological features and help to 
classify different types of manifolds. The kth cohomology group of De Rham of the 
manifold M, H h (M), is the set of closed /c-forms modulo the exact /c-forms. It is a finite 
dimensional vector space, its dimension b^ is called the fcth Betti number. This integer is 
a topological invariant of the manifold M. A beautiful result || [Tj| is that the alternating 



sequence of Betti numbers is the Euler-Poincare characteristics of the manifold, namely 

X(M) = jZ{-l) r b r (23) 



r=0 



This shows that topological properties, as x, can sometimes be calculated using analytical 
tools such as differential forms. 



2.4 Riemannian Manifolds and Absolute Differential Calculus 
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2.4.1 Riemannian Manifolds 



A Riemannian Manifold is a differentiable manifold M with a local scalar product defined 
on each tangent space T P M. If v and w are two tangent vectors at the same point p then 

(v\w) p = v.w = gij(p)v l w 3 . (24) 



In particular the norm of v is given by the square-root of gij(p)v°v- 

The quantities gij(p) are the local components of the metric and are called the metric 
tensor. They allow to compute the length £(7) of a curve 7 parameterized as in (|1|). We 
recall that 4? is the tangent vector to 7 at the point p = j(t). 

V 2 / • -X 1/9 

'i / dx l dx 3 \ ' rW) 




/ .^((-1-),) = / / /f M^^J =X rfs (25) 

The infinitesimal arc length is denoted by ds. In classical books the metric is written as: 

ds 2 = gijdx l dx J (26) 

A geodesic curve between two points p and p' of M is a path of minimal length. Us- 
ing variational calculus, one can find the Euler-Lagrange equations for a geodesic curve 
(xi(s), x n (s)) parameterized by its arc length s (this is a good exercise!). One obtains 
a system of differential equations: 

dv +r} ^ = (27) 



(is 2 ■ ?fe ds ds 
where the Christoffel symbols Tj k are given by 

i.e. the matrix g tl is the inverse of the metric tensor. 



2.4.2 Covariant derivative. Connexion and curvature forms 

We have emphasized that on a general manifold there is no way to compare tangent vectors 
at different points. However, given a smooth vector field, a rather natural question to ask 
is 'What is the infinitesimal variation of a vector field Y at point p if one moves in the 
direction of a vector V tangential at the point p to the manifold ?' A possible way out is 
to immerse the manifold in a larger space where coordinates are defined. But this is not 
intrinsic: it depends on the surrounding space. 

One needs an additional structure, called a connexion, that allows to compare tangent 
vectors at different points and to differentiate vector fields (or more generally tensor fields). 
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This method was invented by Levi-Civita and called absolute differential calculus. It is 
possible to give a general definition of a connexion without refering to the metric, as these 
concepts are independent. However we shall restrict ourselves to Riemannian manifolds. 

A Riemannian manifold can canonically be endowed with a connexion provided by the 
Christoffel symbols T l - k (p8|) . Consider a vector field Y on the manifold. One defines its 

covariant derivative V V Y along the direction of the vector V tangential at the point p to 
the manifold by 

V V Y = + Y) k V^Y k J (29) 

where e, = -£j is the i th vector of the basis (§) of the tangent space T P M at the point p. 
The quantity V V Y is a vector, tangent to the manifold M at the point p. It represents 
the total rate of variation of Y along the direction of V. The first term in ( p9[ ) is nothing 
but the convective derivative of Y along the direction of V. This term is familiar in 
hydrodynamics (e.g. in Euler and Navier-Stokes equations). The second term represents 
the derivative of the vectors e\ of the local basis of T P M along the direction of V. Such 
a term is familiar from mechanics when one uses non-cartesian coordinates (e.g. polar or 
spherical), it reflects that the coordinate system is local and point dependent. One can 
rewrite fl2~9D as follows: 



VyY = V v Y i e i = dY i (V)e i + Y k V v e k 
with V v e k = V ] V ej e k 

and V ej e k = Y) k ei (30) 

It is analogous to the well-known expression in mechanics of the derivative of a vector in 
the rotating frame 

dv d , , . dv => 



*Parallel Transport: The formula (|29"D for the covariant derivative VyY allows to 
compute the variation of a vector field along a given direction. Inversely, from a tangent 
vector Y at a point p and a curve 7 going from p to pi, one can construct from (p9|) a 
vector field Y along the curve 7 such that Y(P ) = Y . Such an operation is called the 
parallel transport of Yq along 7 . The vector field Y is determined by 

= WyY with V = -r . 
v dt 

It is therefore obtained by solving the system of differential equations: 

dxi dt 3k dt dt Jk dt v ' 
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with initial condition Y(0) = Yq. If one considers the special case Y = V = ^ where s is 
the arc-length of the curve 7, equation fl3~T|) becomes identical to the formula defining a 
geodesic (|27|). The equation of a geodesic is thus VyV = with V = In other words 
a geodesic is a curve which is parallel transported along itself, as this is well-known for 
straight lines in Euclidean spaces. 



Remark: We did not give any explanation for the appearance of the Christoffel symbols 
in the formula ([30]) for the covariant derivative of the vectors of the basis. This is in fact 
the heart of Levi-Civita's construction. We encourage the reader to develop an intuition 
on this formula through his/her own readings D], TJ]. A nice (and historical) approach 
is provided by surface theory in R 3 . A tangent vector V at a point p of the surface is 
transported to an infinitesimally close point p 1 by first bringing it parallel to itself from 
p to p' . As V has no reason to be tangent to the surface at p', one projects V on the 
tangent plane at p 1 and obtains a new vector V tangent at p'. The variation of V through 
this operation provides the formula for the covariant derivative (this is explained in ||11|| ). 
Another approach is via 'Tensor Calculus': the first term in the r.h.s. of (|29"D is not a 
tensor, so one has to add to it something to have a tensorial derivative. The Christoffel 
symbols are not tensors themselves but adding them to the first term turns the sum into 
a tensor. This approach is well explained in ||15| , or in books about General Relativity 

H0. 



We now give a more compact expression for the covariant derivative. One first notices 
that V^Y is linear in V, i.e. VF can thus be viewed as a linear operator which acting 
on V produces another tangent vector at p. Equivalently, VY is tangent-vector valued 
1-form and can be written as 

vy = (d Y i + r jk dx ] Y k ) a . 

From this expression we see that V itself acts as a linear operator on the vector Y. 
Recalling that a matrix L acts on Y as LY = L\Y k ei, we write: 

V = d + u (32) 

where uj is a matrix whose coefficients are 1-forms. Precisely u has matrix elements 

o>l = ?) k dxi . (33) 

The matrix- valued 1-form u is called the connexion 1-form and V the covariant deriva- 
tive. The curvature Q of a Riemannian Manifold is defined by the matrix valued 2-form 

Q = did + uj A uj (34) 
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i.e. for two vectors X,Y on the manifold, Q(X, Y) is a matrix. Notice that in general the 
connexion 1-form can be written as 



uj = Ajdx 3 (35) 

where the AjS are matrices. Then, wAw = J2i<j dx l A dx 3 is non zero when the matrices 
do not commute. 

Exercise: Find in a book the formula of the Riemann curvature tensor (see e.g. Spivak 
13 1, Vol II, Chap 4, pl89) in terms of the Christoffel symbols. Using (|33|) , verify that Q 



defined in ( |34| ) is indeed the Riemann tensor. 
2.5 *The Laplacian 

We saw that the spaces of fc-forms and (n — A;)-forms over a manifold of dimension n have 
the same dimension. It is possible to define a duality transformation *, called the Hodge 
star, between these two spaces. In the Euclidean space R 3 the Hodge duality is given by: 

* 1 = dx A dy A dz 

*dx = dy Adz , *dy = dz A dx , *dz = dx A dy, 

*(dxAdy) = dz , *(dy A dz) = dx , *(dzAdx) = dy 

*(dx AdyAdz) = 1 (36) 

More generally the duality * operation on a manifold depends on the local metric, and is 
given by: 

* (dx h A dx i2 A ... A dx ik ) = \g\ 1/2 dx ik+1 A ... A dx in (37) 

where (ii,i2, ik, ik+i, ■■■An) is a permutation of (1,2, ...,n) with positive signature and 
\g\ is the determinant of the metric (gij) (i.e. |g| 1//2 is the volume element). If is a 
fc-form, one can check (exercise!) that * * <\> = (— l) fc(n ~ fc )0 . 

The Hodge duality allows to define a scalar product in the space of k-forms by 



(4>k,i>k) = 4>k a *4> k (38) 

Jm 

where the integral is over the manifold. 

Exercise: Show that the scalar product of two 1-forms = A\dx + B\dy + C\dz and 
■0 = A 2 dx + B 2 dy + C 2 dz in R 3 is 



(0, = / (AxA 2 + B X B 2 + dC 2 )dxdydz. 



The exterior derivative d maps fc-forms onto (fc+l)-forms. One defines the adjoint exterior 
derivative 5 which maps /c-forms onto (k — l)-forms via (<fik,dipk-i) = (^^n) • 
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Exercise: 

1. Show by integrating by part that in a n-dimensional vector space 

S= (-l) nk+n+1 * d* . 

2. Deduce that 5 2 = 0. A co- exact fc-form u satisfies the global condition u = 5<p where 
<p is a (k + l)-form. A co-closed fc-form u satisfies 5uj = . 



The Laplacian is an operation which takes &;-forms onto /c-forms and generalizes the 
usual Laplacian on functions. It is defined by 

A = (d + 5) 2 = d5 + 5d (39) 



Exercise: In R 3 show that if / is a function, one has Af = — (d x 2 f + d y 2 f + d z 2 fj. 

A fc-form ujk such that Au>k = is called harmonic. For a smooth enough £;-form, a 
necessary and sufficient condition for harmonicity is to be closed and co-closed: 

Auk = <^=^> du k = and Suk = . 

Hodge's theorem: On a compact manifold without boundary any fc-form io^ can always 
be decomposed into the sum of an exact form da^-i, a co-exact form 5/3k+i an d a harmonic 
form lk Big: 

uj k = dot k _ x + 5p k+1 + 7 fc (40) 

This very important result is well-known in 3d vector analysis of electromagnetism of 
continuous media as the Helmholtz decomposition according to which for a closed and 
compact manifold or for a controlled growth of the fields at infinity, any smooth vector 
field B can be decomposed as B = V0 iVxM + ff, where H has both vanishing curl 
and divergence. 

Remark: if uo k is closed then d5(3 k+ \ = 0, which implies, as d and 5 are adjoint, that 
SPk+i = 0. Hence u k = da k -i + 7fc . Therefore u k and ^ k belong to the same cohomology 
class. Denoting by Harm fc (M) the space of harmonic k- forms on M, the Hodge theorem 
establishes an isomorphism between the two spaces H k (M) and Harm fc (M) and according 
to (^), the Euler-Poincare characteristic is given by 

n 

X (M) = ^(-l) fe dim Harm fc (M) (41) 

fc=0 

This relation shows that properties of the Laplacian probe the topology of the manifold 

Mm- 
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2.6 Bibliography 



We do not intend to be exhaustive, we just list books we used and give some (subjective) 
comments. A classical and outstanding introduction to Geometry in general is |^TJ. A 
'modern-classic' more inclined toward topology and analysis is ||22|| . A good way to start 
with differential geometry is to study curves and surfaces. A very pleasant book on these 
topics, full of results, figures and historical comments is Some books begin with 

curves and surfaces and then introduce the general concept of n dimensional manifolds. 
This can be a very helpful point of view since it is not too formal [IT, ITS], |23J. More 



recent works are |25| with emphasis on applications, [^] more inclined towards algebra 
and [^UJ devoted to analysis on manifolds and well adapted to theoretical physicists. An 
exhaustive treatment is to be found in H |9], [0J and the following volumes. We owe a lot 
to the review paper [26 written for high energy physics and Yang-Mills theories. For a 
pedagogic introduction to differential forms, see [p7|] . 



3 Fiber bundles and their topology 



3.1 Introduction 

The concepts of connexion and curvature that we have defined for Riemannian manifolds 
can be extended to a more general structure called a fiber bundle. A fiber bundle is a 
manifold X that locally looks like the product of two simpler manifolds. For instance the 
torus T 2 = S 1 x S 1 is globally the product of two 1-dimensional circles S l . A cylinder is 
also a global product S 1 x [—1, 1]. A Mobius strip is locally the product of an arc of S 1 by 
[— 1, 1] but not globally. More precisely, a fiber bundle is a triplet (X, M, 7r) where X is a 
manifold (the total space), M (the base space) a submanifold of X and n (the projection) 
a smooth function from X to M: 



7T : X M 

x i — > n(x) = p (42) 

The inverse image 7t~ l (p) of any point p e M is called the fiber above p and is isomorphic 
to a given manifold F. Intuitively, a fiber bundle is a collection of identical (isomorphic) 
manifolds F which depend on a parameter p belonging to the base manifold M. When 
the fiber F is a vector-space R™ or C n the bundle is called a vector-bundle, (see e.g. f28|j). 



A section s of a fiber bundle is a smooth function which associates to each point 
p G M a element of the fiber above p. Therefore one has tt o s = IcIm where IoIm is the 
Identity function in M. 



Examples: 

1 . The product Mi x M 2 of two manifolds Mi and M 2 is a fiber bundle called the trivial 
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bundle. One can take either M\ or M2 to be the base and the other one to be the fiber. 

2. The total tangent space of M is a vector bundle also called the tangent bundle. The 
base space is the manifold M itself and the fiber above the point p is T p M which is 
homeomorphic to R n . The projection is exactly the one defined in section 2.2.1. And 
what is a section of the tangent bundle? A vector-field! 

3. More generally the space of k- forms is also a vector bundle over the manifold M. 

The fiber bundles described here are assumed to be locally trivial: for any point p of 
the base M there is a neighborhood U p such that the restriction of the bundle over U p , 
i.e. 7r _1 (t/p), is homeomorphic to U p x F . 

3.2 Local symmetries. Connexion and Curvature 

In physical applications one considers vector bundles where the base space is the physical 
space and the fiber is a representation of a continuous symmetry group, i.e. a Lie group 
G with Lie algebra Q. At each point, the order parameter is an element of the fiber over 
that point and a local symmetry group acts upon it. The order parameter is precisely a 
section of a vector bundle. It is natural to ask about the variation of the order parameter 
when one moves from one point of the base to another. However, although all fibers are 
homeomorphic to the same vector space, there is no intrinsic way to identify two fibers 
over two different points p and p' of the base. In other words, given a section s, the two 
vectors s(p) and s(p') can not be compared a priori. This is exactly the same problem 
that we encountered for vectors fields, which are sections of the tangent bundle. The 
extra-structure needed to compute the variation of a section s along a given direction 
tangent to the base is called a fiber bundle connexion and is again denoted by V. The 
expression of V is obtained by reinterpreting the formula fl32| ) : 

V = d + u . (43) 

Here d is the usual operation of differentiation of functions and u is a matrix-valued 
1-form (|3^) on the base space M that can be written 

uj = Aidx 1 

where (x 1 ,...,x n ) is a local system of coordinates on the base M. We impose that u 
represents at each point an infinitesimal transformation of the group G i.e. the matrices 
Ai belong to the Lie algebra Q. The curvature Q of the connexion is also given by the 
same formula as above (|3~H): 

fl = duo + uj A uj = duo + - [u, u\ . 

Important examples: 
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1. For the tangent bundle, we take the symmetry Lie group to be the group of invertible 
matrices GL(n, R). The associated Lie algebra is just the group of all n x n matrices. 
Connexions on a Riemannian manifold are a special case of fiber bundle connexions. 

2. We study now a 2-dimensional physical system M = R 2 where at each point an order 
parameter (or wavefunction) is defined and is a complex number x/j = \ijj\e lx . Hence the 
fiber is F = C. The Lie group for Maxwell electromagnetism is U(l) and the Lie algebra 
is iR. Then, the connexion 1-form (or gauge potential) can be written uo = —iA where 
A = A x dx + Aydy, A x and A y being real functions (lxl matrices). Since the Lie group 
is abelian, the curvature reduces to Q = duo = —iB where B = (d x A y — d y A x )dx A dy is 
nothing but the magnetic field. The covariant derivative given by ( J4"3J ) is then V = d — iA. 
This is indeed the differential operator that appears when one studies the Schrodinger 
equation in a magnetic field, or the Ginzburg-Landau model for a superconductor. 

3. In the Gross-Pitaevskii description for rotating superfluid A He the connexion is the 
Coriolis term Q x r, Q being here the angular velocity. 

4. The Yang-Mills connexion is obtained by the same construction as before with a non- 
abelian Lie group (typically SU(2)). The reader is referred to the original paper of Yang 
and Mills |[29| . The theory of fiber bundles is the right setting for gauge symmetries in 
high energy physics |30], [jjfl . 



3.3 Chern classes. 

The cylinder and the Moebius strip have both S l as base space and [-1,1] as fiber, but 
they are not homeomorphic. One would like to classify different types of fiber bundles 
with given base and fiber. Characteristic classes are cohomology classes of the base 



space that are topological invariants of the fiber bundle |32 . 

We shall study the Chern classes for fiber bundles with Lie group G = GL(n, C), 
the group of invertible complex nxn matrices. The Chern- Weil theorem gives an explicit 
construction of these classes starting from a connexion 1-form u and the associated cur- 
vature Q. This theorem also proves that the invariants obtained do not depend on the 



chosen connexion (see |33| for a most enlightning exposition) . 

We recall that the curvature Q is a 2-form whose coefficients are nxn matrices with 
complex coefficients. A polynomial P(£l) is said invariant if for any matrix g in GL(n, C), 

P(Q) = P(g- l ttg) . 
For example, det(l + j-O) and tr(e^) are invariant polynomials. 

Theorem (Chern- Weil): If P is an invariant polynomial andfl curvature 2-form, then 

(i) PiP) is a closed differential form (dP = 0). 

(ii) If u' is another connexion on the same fiber bundle and Q' the associated curva- 
ture, then there exists a form Q such that P(O') — P(£l) = dQ . 
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Property (i) shows that P(fi) defines a cohomology class of the base space M and 
(ii) proves that this class does not depend on the chosen connexion fi. P(fi) is called a 
characteristic class. In particular, integrals of P(fi) over cycles (i.e. submanifolds without 
boundaries) will provide topological invariants |33], |34] . 

This construction can be generalized to other Lie groups and the polynomial P has a 
different expression according to the associated the Lie algebra. Here, we shall consider 
only the Chern classes defined from 

P(fi) = det(l + ^-fi) = 1 + ci(fi) + c 2 (fi) + ... (44) 

where Cj(fi) is a scalar- valued 2z-form called the i th Chern class. From the expansion of 
the determinant we obtain the expressions 

c = 1 

c i = ^ tr (^) 

c 2 = g-jCtrfi Afi -trfi A trfi) (45) 

and Cj = when 2i is greater than the dimension of M. Since c?P = 0, each Chern form is 
closed as well: dcj(fi) = . The i th Chern class belongs to the cohomology group H 2l M . 
A remarkable fact is that the Chern forms define integer cohomology classes: the integral 
of q(fi) over any 2i-cycle C (i.e. any oriented submanifold C of M, of dimension 2i and 
without boundary) is an integer, depending upon C but independent from the connexion 
1-form uj over the fiber bundle: 

f CiiSl) =neZ (46) 
This integer n is a topological invariant of M called a Chern number. 

Remark: if P(f2) is homogeneous of degree r (see e.g. (|45D ) an explicit formula can be 
given for the form Q that appears in the Chern- Weil theorem: 

Q = r f 1 dtP(u' - Q t , fit) (47) 
Jo 

where fi 4 is the curvature associated to the interpolating connexion = tu' + (1 — t)u. 
Examples: 

1. Consider the U(l) abelian bundle described by the curvature fi = —iB where B is the 
magnetic field. On a two dimensional base space, we have from ([44]) 

p<n) = i + | 
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and the only non-zero Chern class is C\(B) = j-. Then, if M is a closed surface without 
boundary 

— f B = n e Z (48) 
2n Jm 



corresponds to a well-known flux quantization condition |36[ . 

2. More generally, on a manifold M, all the fc-fold exterior products f2,f2Af2,f2Af2Afi 
etc... (k < |dimM) lead to topological invariants. If dimM = 4, only the two first are 
relevant 

Q oc &B.d§ = neZ 

M J 

QAQocf E.Bd 4 x = m e Z (49) 

M J 



To end this section, we discuss the issue of what makes certain physical quantities 
topological numbers and some others not. To that purpose, we consider two examples 
taken from geometry and electromagnetism. 

Example 1: Consider the 2-sphere S 2 of radius R. Its curvature is K = Then the 
area of the sphere J s2 dS = AttR 2 is a metric invariant. But J s2 fi = J S 2 KdS = 4tt is a 
topological invariant and is conserved by smooth deformations of the sphere that change 
the metric. 

Example 2. For electromagnetism, the electric charge inside a closed surface S is given by 

— / E.d§ = — f V.EdV . 

47T JS 47T J 

This statement is equivalent to V.E = p (charge density). The conservation of the electric 
charge in time is 4 § s E.d~§ = 0. It is equivalent to the Maxwell-Faraday equation. These 
two relations do not refer to topological invariance although it is sometimes asserted that 
a topological invariant is obtained by integrating a total divergence. The electric charge 
is a counterexample of this statement. 

3.4 Manifolds with a boundary: Chern- Simons classes. 
3.4.1 The Gauss-Bonnet theorem 

We recall some facts concerning a surface S embedded in R 3 [P5[ . To investigate geo- 
metrical properties, Euler had the idea to study curves drawn on S. More precisely, he 
considered normal sections at a point p, i.e. plane curves obtained by intersecting the sur- 
face S with its normal plane at a point p. Euler proved that such all normal sections have 
a curvature k at p limited by two extremal values K\ and «2 called principal curvatures. 
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Moreover, if the principal curvatures are not equal, there is one curve C\ having curvature 
Ki and another C 2 , perpendicular to the first one, having curvature k 2 . If a normal section 
makes an angle with Ci, its curvature is given by tz = K\ cos 2 (p + k 2 sin 2 . 

Exercise: Prove this theorem of Euler (hint: write the local equation of the surface as 
z = z(x,y) and Taylor-expand to the second order [E3|, O). 



The product of the two principal curvatures Ki and k 2 at a point p is denoted by 
K = KiK 2 - The scalar K is the Gaussian Curvature of the surface S at the point p. 
The Gaussian curvature is conserved under deformations that preserve the metric (e.g. 
bending a surface). It is a metric invariant, that only depends on the intrinsic geometry 
of the surface and not on the surrounding space, although the above definition involves 
this embedding space. Gauss proved this fact by discovering an explicit formula for K in 
terms of the local metric g^- and its derivatives. This is the Theorema Egregium (1827), 



a most 'remarkable theorem' [^, [11], [15| 



Suggestion: Do not try to prove the Theorema Egregium. 

Exercise: Check that the Gaussian curvature of a plane is 0. Prove that the Gaussian 
curvature of the sphere 5* 2 of radius R is K = 1/R 2 . Since S 2 and the plane have different 
Gaussian curvatures it is not possible to make an isometric mapping, even local, between 
them. There is no faithful plane map of the Earth that would respect the distances. 



The Gauss-Bonnet theorem JTTJ for a closed surface S establishes a relation between 



the curvature, which is a metric invariant and the Euler-Poincare characteristic of S: 

lJJ s KdS = X (S). (50) 

The quantity j- J J s KdS where dS is the surface element, is sometimes called cur- 
vatura integra. The Gauss-Bonnet theorem establishes that the curvatura integra is a 
topological invariant of closed surfaces. It has numerous and profound implications. For 
example JRJ, as the reader can check, it implies that no metric of negative curvature can 
be defined on sphere, or that no metric with strictly positive or strictly negative curvature 
can be defined on a torus (for a torus, x — 0). 

One should notice the similarity between the theorem of Gauss-Bonnet (|50|) and the 
Chern numbers (f46D defined as topological invariants obtained by integrating a Chern 
class. The curvature K of a surface can be defined as a characteristic class (called the Euler 
class). So in fact these two expressions are not only similar: they have the same origin. 
Chern developed the theory of characteristic classes while studying higher dimensional 



generalizations of the Gauss-Bonnet theorem [35|. 
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3.4.2 Surfaces with boundary 



We now consider that the surface S has a boundary and is oriented. An orientation means 
that an outward unitary normal vector N can be defined in a coherent manner throughout 
the surface (for instance this is not possible for a Moebius strip). 

Let j(s) be a space curve parameterized by its arc-length s, with ds 2 = dx 2 + dy 2 + dz 2 . 
Suppose that 7 is drawn on the surface S. To each point p G S on the curve, one associates 
a local, orthonormal and direct frame (t, u, N) where t = ^ is a unitary tangent vector 
to the curve 7, u is perpendicular to t and belongs to the tangent plane of S at p and iV 
is the normal vector at p. One has in particular 

u = N x t . 

The curvature vector k of the curve 7 is defined by 

As t is unitary, the curvature vector k has components on u and iV. Its projection on 
N is called the normal curvature of 7, whereas the projection on u defines the geodesic 
curvature k g : 

k g = k.u= (fx ) .iV (52) 



The geodesic curvature k g vanishes for a geodesic curve of the surface S [£3, |15| . We can 
now give the generalization of the Gauss-Bonnet theorem if the the boundary OS is not 
empty. As the curvatura integra is not an integer in general, a boundary contribution 
must be added |TT|: 

hi ' L Kis+ lL k ' dl =* (s) (53) 

where K and k g are respectively the curvature of the manifold S and the geodesic curvature 
of its boundary. 

Example: A spherical cap in the north hemisphere, containing the north pole and limited 
by circle of latitude #0, is topologically equivalent to a disk i.e. to a triangle. Hence its 
Euler-Poincare characteristic is equal to 1. The geodesic curvature of the circle of latitude 
80 is kg = -^cotan^o- The Gauss-Bonnet theorem (|53|) gives indeed 

Area of the cap 1 „ „, . „ .„ 

—— l - + —2nRk g sin 6 = 1 54 

2irR 2 2n 

For 9q — I , i.e. at the equator (half-sphere case), the boundary is a geodesic and k g = 
so that only the first term in (^) contributes. 
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3.4.3 Secondary characteristic classes 

We showed that the Gauss-Bonnet theorem need to be generalized to incorporate the 
contribution of a boundary. Is it possible to modify the characteristic classes theory 
accordingly? The answer is again contained in the Chern-Weil theorem. Since P(Q') — 
P(O) = dQ, the integral of dQ vanishes if the manifold has no boundary. But if there is 
a boundary dM, then by Stokes' theorem 



dQ= Q (55) 

M JdM 

needs not be zero and corresponds precisely the contribution of a 'geodesic curvature'. 
Hence ( |46"D rewrites [34]] : 



f Ci(fi) - / Q(uJ,uJ )=ne Z (56) 

JM JdM 

where uj is the connexion 1-form associated to the curvature Q and ujq is a well chosen 
connexion that compensates boundary effects. The forms Q associated to each Chern 
forms Ci are called Chern-Simons classes or secondary characteristic classes. We now 
present some explicit calculations. The characteristic class associated to Q = duo + uj A uj 
for the Yang-Mills case is tr(f2 A Q). Consider another connexion uj' and the interpolation 
between uj and uj' defined by 

uj t = tu)'+{l-t)u (57) 
Q t = duJt + oj t A uj t 

for t G [0, 1]. Define a = uj — uj' so that by (|47|) , Q is given by 

Q(uj, uj') = 2 f 1 tr(a A QMt (58) 
Jo 

Since Q t = Q — tda + t 2 a A a — ta A uj — tuj A a, then 

Q(uj, uj') = 2tr(a A Q a A da H — a A a A a; a A a A uj a A uj A a) . 

2 3 2 2 

The cyclicity of the trace together with the property ([U|) of the wedge-product provides 

2 

Q(uj, uj') = tr(2a AQ — a A da — 2a A uj A a -\ — a A a A a) . 

3 

Taking uj' = we recover a famous expression for the Chern-Simons connexion [Q namely 

2 

Q{uj) = tr(co> A duj + -uj A uj A uj) (59) 
o 
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When applied to the abelian connexion U(l) on a domain M of the 2d plane, we obtain 
Q(A, A') = A — A' and the curvature is related to the magnetic field Q = —%B. The vector 
potential A 1 to be equal to V%, where \ is the phase of the order parameter ip = \tp\e lx . 
Then, according to (j56f ) we have 

j B+ <f (Vx - A) = n e Z (60) 

JM JdM 

which corresponds to the well-known fluxoid expression. By analogy with the Gauss- 
Bonnet theorem for a surface with boundary (j53|) , we can say that B is a curvature and 
that (Vx - A), which we shall interpret in section 5, plays the role of a geodesic curvature. 



3.5 *The Weitzenbock formula 



The definition of the Laplacian on forms fl39|) is not always adapted to describe problems 



we aim to solve. For a quantum particle moving in a magnetic field described by a U(l)- 
connexion or for the related Ginzburg-Landau equation, the Hamiltonian H (or the free 
energy) is given in terms of the covariant derivative V = d — iA and of its adjoint V* with 
respect to the scalar product on 1-forms defined in (|38|), e.g. H = |V*V in units where 



both the mass of the particle and h are set to one. Thus the question arises to relate the 
covariant derivative to the Laplacian. 

We first denote by D the operator D = d + 8 so that the Laplacian A = dS + Sd = 
D 2 . On a flat space and in the dual basis e$ = dx l , we have Df = J2i e idif so that 
D 2 f = — J2i 9ff coincides with the usual Laplacian on functions. But more generally for 
a Riemannian manifold a curvature term, the local basis vectors are also dynamical 
variables so we write Df = J2i QVj/, where Vj is the covariant derivative along the 
z-direction, then, 

D 2 f = -Y. V?/ + E e;* (ViV, - V.V,)/ 

i j<i 

The second term in the rhs corresponds precisely to the curvature K so that we have the 
general expression 

D 2 f = V*Vf + Kf (61) 

known as the Weitzenbock formula J2^, |37J . 

Coming back to the example of a quantum particle in a magnetic field, we can identify 
the curvature K with the magnetic field. For the case of a constant magnetic field, the 
X-term in (|6"ID is just a constant so that the Hamiltonian coincides with the Laplacian. 
Therefore, on a compact manifold without boundary, the Schrodinger equation is geomet- 
rical and admits topological invariants, namely the total quantized magnetic flux. This 
is not true anymore in the presence of boundaries. For a non-uniform magnetic field, 
the curvature K becomes a local function of the coordinates and the problem defined 
by the Hamiltonian H = |V*V has no geometrical features. To recover a geometrical 
formulation, Aharonov and Casher proved that a Zeeman term (^a.B), playing the 
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role of K in (|6"TD, must be added. Similar non-geometrical features appear for the Pauli 
or Dirac equations. These equations admit zero-modes solutions of the type Df = 
whose number n is an Atiyah-Singer Index. But this number does not correspond to a 
Chern number of the type ( f£B|) since the magnetic field is not a dynamical variable of 
the problem but only a parameter and its magnetic flux is not quantized. To obtain 
an example of a U(l)-connexion with topological invariants, we consider now the case of 
a two-dimensional superconductor described by the Ginzburg-Landau equations, where 
both the order parameter and the magnetic field are dynamical variables. 

4 The dual point of Ginzburg-Landau equations for 
an infinite system 

The existence and stability of vortices in superfluid or superconducting systems have been 
mainly studied for the case of infinite systems or in a limit where boundary effects do 
not play an essential role. Among the large variety of methods available to study vortices 
in superfluids or superconductors, we choose to work in the framework of the Ginzburg- 
Landau expression for the free energy. We consider the case of a finite 2d superconducting 
bounded domain (a billiard) and study the existence and stability of vortices. The su- 
perconducting state is characterized by a complex order parameter. For infinite systems, 
nonlinear functionals given by Ginzburg-Landau, Gross-Pitaevskii or Higgs expressions 
admit vortex like solutions. These solutions are characterized by topological numbers e.g. 
the number n of vortices. How can these results be extended to finite size systems? Is 
there a mechanism by which boundary conditions may allow to select a state with a given 
number of vortices ? 

4.1 The Ginzburg-Landau equations 

The Ginzburg-Landau equations describe a superconducting billiard if both the order 
parameter and the vector potential have a slow spatial variation. The expression of the 
Ginzburg-Landau energy density a is 

2e - B 2 

a = a + a 2 \ip\ 2 + a 4 |^| 4 + a^V - i— A)^\ 2 + — (62) 

TlC 87T 

where ip = \ip\e lx is the complex- valued order parameter, B is the magnetic field and the 
Oj's are real parameters. Defining |$9| £ 2 = j^y, A 2 = ^v^2(ff ) 2 q 2 | ' ^ e dimensionless 
free energy T is 

T= l \\B\ 2 + k 2 \\ - H 2 | 2 + |(V - %A)^\ 2 (63) 
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where ijj is measured in units of i/jq = y j^j-, B in units of 4^2, and the lengths in units 
of A\/2. The numerical factor a/2 is for further convenience. The ratio k = ^ is the only 
free parameter in (|63|) and it determines, in the limit of an infinite system, whether the 
sample is a type-I or type-II superconductor p9 |. The integral is over the two-dimensional 



domain f2 of the superconducting sample. The Ginzburg-Landau equations for the order 
parameter if> and for the magnetic field B = Vxi are obtained from a variation of T . 
They are nonlinear second order differential equations. Their solutions are not known 
except for some particular cases. 

4.2 The Bogomol'nyi identities 

For the special value k = the equations for ip and A can be reduced to first order 
differential equations. This special point was first used by Sarma in his discussion 



of type-I vs type-II superconductors and then identified by Bogomol'nyi |30| in the more 
general context of stability and integrability of classical solutions of some quantum field 
theories. This special point is also called a duality point. We first review some properties 
of the Ginzburg-Landau free energy at the duality point. We use the following identity 
true for two dimensional systems 

\{W-iA)i)\ 2 = \V^\ 2 + V x j + B^f (64) 

where j = Im(ip*'Vil ! ) — \4>\ 2 A is the current density and the operator T> is defined as 
D — d x + id y — i(A x + iA y ). This relation is a relative of the Weitzenbock formula flBTP . 
At the duality point k = the expression (|63| ) for T can be rewritten using (p4]) as 

F= ^Qlfi-l + l^lT + l^l 2 ) +fjj + A).dl (65) 

where the last integral over the boundary dQ of the system results from Stokes theorem. 
For an infinite system, we impose |^0[ the usual conditions for a superconductor, namely 



|-0| — > 1 and at infinity. The boundary term in ( |55] ) then becomes 

/ (j + A).dl=i {J-+A).dl (66) 
Jan Jan \ip\ z 

This last integral is the fluxoid. It is quantized and is equal to 

/ Vx-dl = 2nn . 

JdCl 

The integer n is the winding number of the order parameter i/j and as such is a topological 
characteristic of the system. Using we see that n is also the total magnetic flux 

through the system: 

' B = n . (67) 
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As we interpreted B as a first Chern class, this relation is similar to (|46|). The extremal 
values of J 7 , are obtained when the two Bogomol'nyi |^0[ equations are satisfied 



Vip = 
B = 

In this case, the free energy T is given by 







2tt 



n 



(68) 



(69) 



Therefore the free energy itself is a topological invariant. The two Bogomoln'yi equations 
can be decoupled and is a solution of the second order nonlinear equation 

V 2 /n|Vf = 2(|Vf - 1) (70) 

which is related to the Liouville equation. 

It should be noticed that the set of equations ( p5| , |7D| ) has been obtained without any 
assumption on the nature of the magnetic field and appears in various other situations, e.g. 



Higgs [ 43) , Yang-Mills |EJ and Chern-Simons [ 44] field theories. It was proven that these 
equations admit families of vortex solutions [4"5|. For infinite systems, it can be shown 
that each vortex carries one flux quantum and that the winding number n is equal to the 
number of vortices in the system. However for an infinite system there is no mechanism 
to select the value of n. It will be precisely the role of the boundary of a finite system 
to introduce such a selection mechanism and to determine n, according to the applied 
magnetic field. 



5 The superconducting billiard 

From now on, we shall study finite size systems in an external magnetic field B e and 
we denote by <p e = nR 2 B e /<j) its dimensionless magnetic flux through the system. The 
question then arises to know if they can sustain stable vortex solutions and what their 
behaviour is, as a function of the applied field. An interesting approach has been developed 
by Bethuel and coworkers |^7|]. They considered the case of a billiard Q without applied 



magnetic field but with the boundary condition for the order parameter = g(6) with 
g(6) = and a prescribed winding number ^ Jff n %d@ = Ji . In the London limit, i.e. 
k — > oo, is 1 almost everywhere but because of the degree n on the boundary, must 
vanish n times in the bulk therefore leading to vortices. An extension of this approach to 
the case where a magnetic field is applied on the system has been proposed in where 
it is shown by a variational argument that vortex solutions have a lower energy when the 
magnetic field is increased. By the same method, it is also possible to discuss the type 
of vortices and their distribution as a function of the geometry of the billiard. Numerical 
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simulations |35j of the Ginzburg-Landau equations for a long parallelepiped in a uniform 
magnetic field show that the physical picture derived for k = namely the existence of 
stationary vortex solutions whose number depends on the magnetic field, remains valid 
for quite a large range of values of n, and the corresponding change of free energy is small 
46f . We shall therefore study the case k = i.e. the duality point and extend the 



previous approach to a system with finite size where boundary effects are important. 
5.1 The zero current line 

In a finite system, there are in general non-zero edge currents and the order parameter is 
not equal to 1 on the boundary. Hence, the identification of the boundary integral in fl65|) 
with the fluxoid (|66j) is not possible anymore, and the free energy can not be minimized 
just by imposing Bogomol'nyi equations fl68|). However, the currents on the boundary of 
the system screen the external magnetic field and therefore produce a magnetic moment (a 
circulation) opposite to the direction of the field, whereas vortices in the bulk of the system 
produce a magnetic moment along the direction of the applied field. Hence currents in the 
bulk circulate in a direction opposite to those at the boundary. If one assumes cylindrical 
symmetry, J has only an azimuthal component, with opposite signs in the bulk and on 
the edge of the system (the radial component is zero since J is divergence free). Thus, 
there exists a circle T on which / vanishes. This allows us to separate the domain Q into 
two concentric subdomains Q = Q\ U Q2 such that the boundary dQ\ is the curve T. On 
dfl\, the current density /is zero, therefore 

an lJl = L W Jl = °" (71) 

Thus one deduces as above that Bogomol'nyi and Liouville equations are valid in the 
finite domain Qi as in the case of the infinite plane. The magnetic flux <&(f2i) is calculated 
using the fluxoid and ([fl]) so that 

*(n 1 )=n-/ f^ = n. 

As before n is the winding number, i.e. § dn Vx-dl = 27m, as well as the number of 



vortices [^9| in Vt\. The free energy in Vt\ is 

= 27m. (72) 

The contribution of Q2 to the free energy is given by (2) and can be expressed using 
the phase and the modulus of the order parameter ip 

m) = / Wl) 2 + |Vf|Vx - A\* + f + (1 ~f |2)2 (73) 
JQ.2 2 2 
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The boundary conditions for both the magnetic field B(R) and the vector potential 
A(R) adapted to a flat disk geometry are provided by the condition (f> = <f) e . It implies 
that at the boundary B(R) is larger than the applied field B e due to the distorsion of the 
flux lines near the edge of the system. 

5.2 A selection mechanism and topological phase transitions 

The analysis presented in E9| leads for the total Gibbs potential of the billiard to the 



expression 

-^-G(n,(p e ) = ^-F{n,(j) e ) - ^-0 e 2 

Z7T Z7T R z 

= n+^{n- e ) 2 - \(^?{n - e ) 4 - ^0 e 2 (74) 

This relation consists in a set of quartic functions indexed by the integer n. The minimum 
of the Gibbs potential is the envelop curve defined by the equation §^|<£ e = (using 
obvious notations), i.e. the system chooses its winding number n in order to minimize Q. 
This provides a relation between the number n of vortices in the system and the applied 
magnetic field e . In the limit of a large enough y, the quartic term is negligible and the 
Gibbs potential reduces to a set of parabolas. The winding number n is then given by 

» = w.-^ + 5i (75) 

The magnetization of the system M = —§r- is given by 

-M=^(0 e -n)-^0 e (76) 

For (p e smaller that 2 ^ x , we have n = and (— M) increases linearly with the external 
flux. This corresponds to the London regime before the first vortex enters the system. The 
field B\ at which the first vortex enters the system corresponds to Q(n — 0) = Q{n — 1), 
i.e. to B\ = 2iT ^ RX + 2^r2- T ne subsequent vortices enter one by one for each crossing 
Q{n + 1) = Q{ri)\ this happens periodically in the applied field, with a period equal to 
AH = -^2-. This gives rise to a discontinuity of the magnetization AM = 2 "^~\ 



5.3 A geometrical expression of the Gibbs potential for finite 

systems 

For infinite system, with the boundary conditions |^| — > 1 and at infinity, we showed 
that the free energy at the dual point k = 4^ is a topological invariant proportional ( |69D 
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to the fluxoid n, which represents also the number of vortices in the system: 

l?=jB = n. (77) 

This relation is a property of the dual point at which the Ginzburg-Landau functional 
has a geometrical interpretation and admits the quantized magnetic flux as a topological 
invariant. For a finite system the fluxoid quantification can be expressed as 

ii Sj3+ Lw d=n - (78) 

This relation is analogous to the Gauss-Bonnet theorem for a surface with a boundary 



(|53|) or more generally to a topological invariant obtained by summing a Chern class in 



the bulk and a Chern-Simons class on the boundary (pq) . Recalling (pOD, we see that 
the magnetic field B plays the role of a curvature K, and the quantity j^p- = Vx — A is 
similar to a geodesic curvature k g . For a system with cylindrical symmetry, one can show 
p9| that T^p- — n — 4> e . Therefore equation flT^]) can be rewritten as 



Comparing ( [79] ) with ( |77D we conclude that the boundary correction is a functional of 
the geodesic curvature. In the preceeding section, we obtained an explicit formula for the 
function r\ as an even fourth order polynomial in the geodesic curvature. 



This geometric interpretation makes us believe that an expression such as (|79|) is fairly 
general. It could be well suited, as an Ansatz, to describe finite systems which are known 
to have a topological description in the infinite limit; for example, a suitable generalization 
of (|7^) to SU(2) symmetry could describe superfluid 3 Helium in a bounded domain. 

We have presented a geometrical formulation for the Ginzburg-Landau problem which 
we now briefly summarize. 

i) For a certain theory, like Ginzburg-Landau or other functionals of this type, it may 
appear stable singular solutions (e.g. vortices) whose nature is determined by the related 
homotopy groups of the Toulouse-Kleman [|TJ approach. To these solutions are associated 
topological numbers in the Bogomol'nyi limit. 

ii) The existence of topological numbers signals the occurence of a geometrical de- 
scription of the problem. 

What is it good for? First, we notice that topological quantities describe global fea- 
tures of the problem i.e. behaviour in the large by opposition to the local behaviour 
obtained from solutions of differential equations. Then, by identifying physical quantities 
in terms of global topological invariants, we do not need to solve the local equations to 
obtain the behaviour of the system. As a result, we may say that an important goal of 
a geometrical description is to obtain physical quantities in terms of global topological 
expressions. When it is possible, it is very much rewarding. 
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